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Abstract 

We consider a Dirac field in 2 + 1 dimensions with a domain wall like 
defect in its mass, minimally coupled to a dynamical Abelian vector 
field. The mass of the fermionic field is assumed to have just one linear 
domain wall, which is externally fixed and unaffected by the dynamics. 
We show that, under some general conditions on the parameters, the 
localized zero modes predicted by the Callan and Harvey mechanism 
are stable under the electromagnetic interaction of the fermions. 
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1 Introduction 



It is a well-known fact that, in an odd dimensional spacetime, a domain wall 
defect in the mass term of a Dirac field induces a fermionic zero mode local- 
ized on the defect [p]]. This effect is known to occur even in the presence of an 
external gauge field, if the corresponding electromagnetic field is contained 
in the defect hyperplane. Different aspects of this kind of system have been 
studied both for static @, |3| , and dynamical H defects. As far as we know, 
however, possible effects due to interactions between the fermions have not 
been considered for this system. In this article, we shall study the stability of 
this kind of configuration when the electromagnetic interaction between the 
fermions is turned on. That the localization phenomenon should survive this 
interaction is not apriori evident. For example, for a static configuration, 
the Coulomb repulsion between the localized charges could be so important 
as to spread the charge density out over a large region, since the charge den- 
sity due to the zero mode shall induce an electromagnetic field normal to 
the defect hypersurface. On the other hand, we note that our study may be 
thought of as a domain-wall analog of the consideration of the self consistent 
vacuum currents in the presence of vortices [[J . 

This paper is organized as follows: in section [2], we introduce the model 
and derive a self-consistent equation based in some approximations. This 
equation is solved for two different mass profiles in section |^. Finally, in 
section |] we discuss the effects of the non-zero modes and present our con- 



The Euclidean action S, for the system we shall consider, is given by 



elusions. 



2 The model 



s 



Sf + Sg 



(1) 



where 




x ij){x) [ft + it 4(x) + M(x)]if)(x) 



(2) 



is the fermionic action, and 




(3) 
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the Maxwell action, which defines the gauge field dynamics, x = (x ,Xi,x 2 ) 
denote the Euclidean coordinates, and the Hermitian 7 matrices are assumed 
to be in an irreducible 2x2 representation of the Dirac algebra, verifying the 
anticommutation relations {7^,7^} = 2 5^. The complete Green's functions 
can be derived from the generating functional 



+ / d^xlj^A^x) + r)(x)ip{x) + ip(x)r)(x)} 



(4) 

where we included source terms for the gauge and fermionic fields. The 
fermion mass is regarded as an external classical 'field', dependent on the x-i 
coordinate only. We also fix the number of defects to one, by requiring M{x) 
to cross once, at x 2 = 0, say. 

By applying the property that the functional integral of a (functional) 
derivative vanishes to equation (f|), we derive the 'quantum equations of 
motion' 







VA^VtpVip 



5S 



exp{-S[ij,ip]A^ 



+ / d 6 x\j ll {x)A ll {x) +t}(x)i/)(x) + ip(x)r}(x))} 



for Ap, and 



(5) 



VAJ>iJiThf) 



5S 
5ip(x) 



i](x) 



exp{-£[<0,-0; A h 



+ / d 3 x\jy l (x)A /J ,(x) +r}(x)ijj(x) + ip(x)r)(x)]} , (6) 

for ip (the adjoint equation is trivially obtained). Taking the functional 
derivative with respect to rj(y) in @, and putting all the external sources 
equal to zero afterwards, we find that equations (El) and (EH reduce to: 



JJx) 



and 



[^ + ieA(x) + M(x)]ip(x)ip(y)) = 5(x - y) 



(7) 
(8) 
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where: 

J v (x) = ie(^(x)j v i>(x)) (9) 

and 

= dpAv-dyAy. , A = (A^) . (10) 

Equation ([?[) is an inhomogeneous 'classical' Maxwell equation, with the 
average gauge field A^ = (A^) playing the role of the classical gauge field, 
and the average (vacuum) fermionic current as its source. Equation @ 
involves the expectation values (ipip) and (Aipijj). Of course, an exact treat- 
ment would require the use of an infinite set of coupled equations involving 
all the different Green's functions of the system. In order to find a simpler 
and closed system of equations, we make the following approximation: 

(A^xMxMy)) ~ (A^x))ty{x)$(y)) = A,{x) S A (x,y) , (11) 

where we introduced S A (x,y), which denotes the fermionic propagator in 
the presence of an 'external field' A(x), which corresponds to the average 
gauge field. This amounts to a sort of mean field approximation, where the 
gauge field is treated classically. To make the approximation involved more 
explicit, we note that the (exact) three point function appearing in (|TTJ) can 
be written in the equivalent form: 

(A„(x)^(x)^(y)) = JvAA^ix^ + iefi + MY^e-^-^^ (12) 

where 

T F [A] = -logdet[^ + ie A + M] . (13) 

The approximation ([□]) is obtained from ([L2"|) by replacing A by its saddle 
point value. Namely, the approximation amounts to using the (leading) 
saddle point approximation, where the 'action' which is minimized at the 
saddle point is the bare Maxwell action plus an effective contribution IV [A] 
coming from the fermionic determinant. 

Equation ( |TTD is sufficient to close the system of equations, since then (^) 
becomes: 

[ft + ieA{x) +M(x)}S A = 5{x-y). (14) 

It is now important to realize that the average current can be expressed as a 
functional of A, as follows: 



nfp(if>(x)ij){x)) = -ietr[y lt SA{x,x)] . (15) 
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Equation (|7|), together with (|H>D, define a closed system of equations, which 
allows us to find the average gauge field A, and then the current density 
induced in that background. The equation that determines A is obtained by 
replacing by its expression (|T5|) in (^): 

d^F^ u (x) = -ietT[jfj,S A (x,x)} , (16) 

which, in general, and depending on the approximation used to evaluate S^, 
will be a non-linear integro-differential equation. The non-linearity comes 
from the fermionic propagator S4, which is defined as : 

S a p(x,y) = (x,a I V~ l \ y,@), (17) 

where V = (0 + ie A + M). 

We shall now look for particular solutions of the coupled set of equations, 
under some restrictions and simplifying approximations. We shall restrict 
ourselves to static, purely electric solutions, with no electric current (hence, 
no magnetic field). In the Coulomb gauge, the only remaining component 
for the (average) gauge field is Ao, which is determined by the equation 

VV = -ietr[ lo S v (x,x)} , (18) 

where V — Ao. 

Our approach to solve the system of equations shall be to first evaluate 
the fermionic propagator in the external potential V. Then, we shall find 
the corresponding vacuum charge density as a functional of V, and insert it 
into the Gauss law (|18|) to determine V. The resulting V can then be used 
to fix the precise form of the charge density. We will be able to say that 
there are localized modes if the system admits solutions where the charge 
density is confined to a small region around the defect. Of course, we shall 
have to make some assumptions also on the allowed boundary conditions for 
the fields. The choice of these conditions is also part of the kind of ansatz 
used, and also on the amount of generality one wants to introduce into the 
treatment. 

To find the fermion propagator in the presence of the external field V, 
we shall use the perturbative expansion of T>~ 1 in powers of V, namely, we 
decompose V as follows: 

T> = V + V, (19) 

where 

V o =0 + M(x) (20) 
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and 

V = ie-foV(x) . (21) 
Thus, T>~ x is naturally expanded as: 

V v = V^ 1 - V^VVq 1 + V^VV^VVq 1 - ... (22) 

We note that the 'free' propagator Vq 1 includes the mass field and its space 
dependence exactly. This must be so, since the defect changes the spectrum 
of the Dirac field, an effect that cannot be described perturbatively. To find 
the inverse of V , we use the equivalent expression: 

V, 1 = (VlVoy'vl (23) 

which requires finding the inverse of the Hermitian operator 

Ho = VlV . (24) 

This is a much simpler task than inverting V , and it allows one to dimen- 
sionally reduce the problem. To see this, we follow the procedure of @, of 
which we give a lightning review here. First we write: 

V = (a + @)P L + (a^ + })P R , (25) 

where = 7o<% + 71 d\. We define the operators a* and a, that act on 
functions of the X2 coordinate as 

a = d 2 + M a) = —d 2 + M, (26) 

and the projectors Pl, Pr- 

These projectors behave like chirality projectors from the point of view of 
the 1 + 1 dimensional theory which describes the chiral zero mode. This 
decomposition makes it possible to disentangle the dynamics corresponding 
to the £2 coordinate from the coordinates ). The 'dimensional 

reduction' can be seen to arise at the level of the operator TL$: 

H = {h- f)P L + {h- f)P R , (28) 
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where 

h = a) a h = a a) . (29) 

To expand the fermionic fields, we define n and n , eigenstates of the oper- 
ators h and h, respectively. We denote by A^ their (common) eigenvalues: 

h(p n = \ 2 n (fr n , h(j) n = \ 2 n (j) n , (30) 

(<f>n | <Pm) = 8nm, (4>n \ 4>m) = S nm , (31) 

since the spectra coincide, except for A n = 0, and the eigenvalues are of 
course positive. The X n = eigenvalue will, by assumption, be present only 
for h. This will depend of course on the mass profile near the defect, i.e. 
the zero of the mass. Since the sign of A n is arbitrary, we take it positive by 
convention. 

Thus, the fermionic fields can be expanded as: 

^(x,x 2 ) = YM>n(x 2 )^\x) + Mx 2 )^\x)}, (32) 

n 

ftx,x 2 ) = £[^ n) (x)4(x 2 ) +^t\x)ft n {x 2 )}. (33) 

n 

The spinors that carry the dependence on x are defined by: 

<k*) = PL,R^ {n \x) , ^ R {x) = ^ in \x)P R , L , (34) 

where V'iTii denotes a general bidimensional fermionic field (one for each value 
of the index n). In terms of this expansion, the fermionic action becomes: 

s = + E s(n) ( 35 ) 



where S L denotes the action for a chiral left-handed fermion in 1 + 1 dimen 



(0) 
L 

sions, while is a massive Dirac action, also in 1 + 1 dimensions, with a 
mass equal to A n (the sign of the mass is irrelevant in 1 + 1 dimensions). 
Since 

Wo 1 = (h- fy l P L + (h- fy l P R , (36) 
the free propagator becomes 

V, 1 = {h- f)- x P L (ct -p) + (h- fr^Ria - }) . (37) 
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Translation invariance along the x and X\ coordinates suggests the use of 
a potential depending only on i 2 , V = V(x 2 ). To find the propagator in 
configuration space, we need to evaluate the following expression: 

S a p(x, y) = {V^U{x, y) - (V^VV^^x, y) + ... (38) 

with V = V(x2). In the perturbative expansion for the propagator, we insert 
expansions of the identity constructed with intermediate states correspond- 
ing to eigenstates of the operator TCo- Using the fact that each eigenvalue 
A n corresponds to the effective mass of a two dimensional mode, and that 
the lowest mode is massless (the zero mode), it is natural to keep only the 
zero mode in the intermediate states as a first approximation. Note that the 
mass A n of the non zero modes is separated from the zero mode by a finite 
gap whose magnitude is controlled by the profile of the mass near the defect 
(see ref. @). With this in mind, we shall first use the leading approxima- 
tion of keeping just the zero mode, and then make a quantitative evaluation 
of the error involved in this procedure, by including the correction corre- 
sponding to the lowest massive mode. On the other hand, we shall keep the 
full dependence in the potential, namely, we shall use no truncation for the 
perturbative series in V. To implement this approximation, we introduce 
projectors Pq along the zero mode. They are explicitly given by 

Po = h4E^L ) - (39) 

n 

Taking this into account, after some algebra one can show that, in this ap- 
proximation, the propagator is given by: 

Sap(x,v) - <Po( x 2)4>o{y2)(x ,x 1 ,a | - J~ 16 1oV ®? Pl o I yo,yi,P) ■ (40) 

(d + le 7o^o,o) 2 

In this expression there appears the average of V in the zero mode which is 
denoted by: 

V , = (0o|U|0 o ) • (41) 

It is worth noting that this result is approximate in the sense that only 
the zero mode has been included, but all the powers of have been added, 
as it is evident from the non-linear dependence of the propagator on A^. 
The charge density is evaluated by multiplying by 7 , taking the Dirac trace, 
and finally calculating the coincidence limit x — > y. Inserting the result so 
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obtained for the charge density as a functional of the potential into (IB) 
yields: 

-V(x 2 ) = (j)o{ x 2m{x2) / - , . ] — — 77^ • (42) 



dx\ J 2tt (-i k + ie-f V 0fi ) 2 

The momentum space integral has both linear and logarithmic divergences. 
Using a symmetric limit kind of regularization, we see that (|42|) can be ex- 
pressed as: 

d \v(x 2 ) = o (x 2 )4( X2 )^M . (43) 



dx\ 

It is remarkable that, as a consequence of the fact that we are only keeping 
the zero mode, the expression for the charge density becomes linear in the 
potential. This happens in spite of the fact that we have kept all the powers 
of the potential in the fermionic propagator, since the result is a consequence 
of the fact that massless two dimensional QED is exactly solvable ||, with 
the exact fermionic determinant being quadratic in the gauge field. 

We have obtained an integro-differential equation involving derivatives of 
V and its average on the lowest energy mode. To solve it self-consistently, 
we first derive from (f43|) (by integration) an equation for V, depending also 
the average of the potential. Then, as a second step, we shall insert this 
average into (f|^) in order to obtain the explicit profile of the potential as a 
function of x 2 . At this point, it is clear that the existence of a self-consistent 
solution depends on the particular form of the zero modes appearing in eq. 
fl4"3]). This differential equation will have a solution only if the charge density 
is localized in such a way that the integrals involved are well defined. In 
particular, the zero modes need to be localized around the defect. It was 
shown in reference [|IJ that in 2n + 1 dimensions the zero mode has the form 

r]e -J: 2 dyM(y)^ (M) 

where rj is an spinor independent of x 2 - 



3 Examples 

In what follows we will discuss the 
different kinds of mass profiles. 

• Step-like defect. 



possible solutions of eq. (|43|) for two 
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Given a mass of the form: 

M(x 2 ) = A(2Q(x 2 ) - 1) . (45) 

where A is a constant with the dimensions of a mass, and G is the Heav- 
iside function, there is only one zero mode [|2|], which can be explicitly 
written as: 

o (x 2 ) = A^ e ~ A NI . (46) 
In this case, the differential equation becomes 

|k(* 2 ) = iAeVo, e- 2A l-', (47) 

and integrating it twice we obtain for the potential V 

V(x 2 ) = a+^ K e 2 V 0fi e- 2A ^ , (48) 

where a is a constant, to be related later to the chemical potential. 
In this expression we have not included a term that corresponds to a 
constant electric field in the x 2 direction, because it could be eliminated 
by choosing appropriate boundary conditions (such as vanishing density 
of charges at infinity) . 

In order to find a self-consistent solution for the potential we evaluate 
the expectation value of V, which is expressed by fl48|), in the zero mode 



1 POO 

V ,o = a + -eV„, / dx\Mx)\ 2 e" 2 ^ . (49) 
oA J -oo 

Thus, Vo,o is easily seen to be given by 

V ,o = — V- • (50) 

1 16A 

Therefore, the potential written in terms of the zero modes results 

2ae 2 

V(x 2 ) = \Mx)\ 2 16A2 _ Ae2 - (51) 

Notice that the solution is only stable if the electromagnetic coupling 
constant and the mass coupling constant satisfy the bound: e 2 < 16A, 
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which means that the strength of the interaction (repulsion) between 
the electrons cannot be larger than the scale given by the height of the 
defect. We note that 'stability' refers here to the property of having a 
confining potential. We see that in this case, i.e., for an step-like mass 
and keeping only the zero energy mode, there exist a self-consistent 
solution for the fermionic interaction potential. In other words, even 
in the case of interacting electrons, the fermions are localized in the X2 
direction and can only move along the defect. 

The interpretation of a as a chemical potential proceeds from the fact 
that the Gauss law (47), combined with (0), means that the charge 
density of the configuration is 

PM = a e- 2A '^', (52) 

A 1 16a) 

and (by integrating over x%) one sees that the total charge is propor- 
tional to the constant a. 

Linear defect. 

Assuming than the mass can be expanded as a power series in x 2 > for 
small enough X2 we only keep the first order term: 

M(x 2 ) = M'(0) x 2 , (53) 

where we assume M'(0) 7^ being M' the first derivative of the mass. 
For this mass profile we can still find the zero mode by defining [|| 

h = -d\ - M' + M 2 x\ . (54) 

which is an harmonic oscillator Hamiltonian. The lowest energy mode 
is: 

,\M'\ u/4 _Wll x 2 
<t>o{x2) = ( J -) ' e 2 *2. (55) 

7T 

Following the same steps as in the previous example we find that the 
potential can be written in terms of the zero mode as 

V(x 2 )=a + Hdy [ y d z \<t> ( z )\*(-^—), (56) 
jb ja I — e z G 

Wllere roc rx 2 r y 

C= dx 2 \Mx2)\ 2 dy dz\Mz)\ 2 - (57) 

J-00 JB JA 
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Thus we see that also in this case there exists a self-consistent solution 
for the Gauss law, for a charge density localized around the defect. 
However there is a necessary condition for the existence of this localized 
mode. The wave function of the zero mode has to vanish rapidly outside 
the region of the space where the mass can be approximated linearly. 
A quantitative criterion for the validity of this condition can be found 
in reference 0. 

In summary, up to know we have shown the existence of localized so- 
lutions if we keep only the lowest energy modes in the expansion of the 
fermionic propagator. This solution depends on the mass profile, and it is 
non-perturbative in the electromagnetic interaction between the fermions. 
We have neglected the (more energetic) massive modes based on the fact 
that the terms on the action that come from these modes go as tt, where X n 
is the mass of the mode 0. Therefore, for a large and steep enough mass, 
our approximation will be valid. In particular, for a linear defect, the mass 
of the modes is proportional to the slope of the mass profile. Therefore, by 
changing this slope we could make A n arbitrarily large. 



4 Effect of the massive modes 

We shall now study the problem of including one massive mode in our calcu- 
lation in order to check whether there still exist localized solutions or not. It 
will provide also a quantitative idea about the error involved in considering 
only the lowest energy modes. 

We proceed as follows: in the perturbative expansion for the fermion 
propagator (|22"|), we consider only the projection of the operators Vq 1 and V 
onto the two lowest energy modes. Then each factor Vq 1 contributes with 

Z>o 1 ~ 0kMAiPl - 4>i4$Pr + M%Pr ~ M\PPl - MIPPl- (58) 
For an even V(x2), selection rules imply the vanishing of the matrix elements 

^71,771) ^n,rh and Vfi,fh- 

Replacing ( |58|) into ( |22[), and keeping only the non- vanishing matrix 
elements of V, the first order correction in V to the fermion propagator (i.e., 
correction to ( |58|)) is: 
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(Af - foffl ^ + 7^jr [ ^ 7 v °'°^ Pl - (59) 

Notice that in this case it is not possible to obtain a non-perturbative ex- 
pression for the fermion propagator due to the fact that we are taking into 
account massive modes as well as the massless one. In order to write the 
Gauss law we need to compute 

+ ^^ [(Mg - S^oa] + ^[(286 - &)VU (60) 

Taking the Fourier transform in the above expression and regularizing the 
integrals by a symmetric limit, the Gauss law becomes 

d 2 e 2 V\ i e 2 V- - 

V(x 2 ) = <pi{x 2 )4>\(x 2 ) — — ^ + 0i(x 2 )0i(x 2 ) — (61) 



Thus we have obtained a differential equation whose solution will depend on 
the localization properties of the fermionic modes around the defect. 

In the case of a mass that can be approximated by a linear function of x 2 
near the defect, it is simple to check that @ 

4> n = (}) n+ i, (62) 

therefore there is only one zero mode, and the Gauss equation becomes 

JL-V(x 2 ) = <pQ(x 2 )<t>l(x 2 ) e + <pi(x 2 )<p\(x 2 ) e (63) 
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Integrating this expression we find 

V(x 2 )=a + (^) ["dy [" dz\Mz)\ 2 + (^) H dy P dz\<h(z)\ 2 - 
2 J b J A 2 Jb J a 

(64) 

Once again, we look for the self-consistent solutions for the expectation values 
of the potential. When computed on the two lowest energy modes, they are 
given by the solution of the equations: 

V iti = a + (^)Ao + (^)Ai, (65) 
where i = 0,1 and are: 

/oo rX2 ry 

dx 2 \(t )i {x 2 )\ 2 / dy / dz\^(z)\ 2 . (66) 
-oo JB J A 

Solving (|65l ) we obtain 

t/ - 1 - e 2 D 10 + e 2 D m 

V °'° ~ * (I - eW 01 )(l - eW w ) - eW llDoo > (67) 



Via 



1 - e 2 An + e 2 D u 



(1 - e 2 An)(l - e 2 D 10 ) - e*D u D 00 
We have found that, in the case of a linear mass, there exist a self-consistent 
solution of the Gauss equation to first order in the interaction potential, if we 
include apart from the zero mode, one massive mode. Notice that, for a linear 
mass around the defect, <p n and cj) n are harmonic oscillator eigenstates. Far 
enough from the defect, the eigenstates decay exponentially (as a Gaussian 
function), ensuring that the charge density is localized around the defect in 
such a way that there is a solution for the Gauss equation. Obviously all the 
caveats regarding the range of validity of approximating the mass by a linear 
function, that we mention in the previous case, must be taken into account 
here. 

Summarizing, we have considered a Dirac field in 2 + 1 dimensions with a 
domain wall like defect in its mass, minimally coupled to a dynamical Abelian 
vector field. The mass of the fermionic field is assumed to have just one linear 
domain wall, externally fixed and unaffected by the dynamics. In the absence 
of electromagnetic interactions among the fermions, it is a well known fact 
that localized zero modes exist on the defect 0]. We have studied here the 
effect of the fermionic interactions on these modes showing that, under some 
general conditions on the parameters, the localized zero modes stable under 
the electromagnetic interactions of the fermions. 
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